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1 Introduction 



In a recent paper Kong and Ravndal examined the problem of the decay 
width of the pionium atom — 7r + 7r~ — using methods of effective field theory. Q 
Specifically they employed a quasi-local approximation to the pi-pi interac- 
tion and evaluated the width using simple second order perturbation theory 
to obtain the canonical answer 

r^-^h^l^o)! 2 (i) 

in terms of the 7r + 7r~ — > 7r°7r° scattering length a 12 . Here k\ = 
y / 2m 7r o(m 7r + — m^o) = y/2m 7r oAm 7T is the center of mass momentum in the 
final state — 7r°7r° — system and |\I/(0)| 2 = m^ + a 3 /87r is the wavefunction at 
the origin. In a subsequent paper they examined the problem of pp scattering 
within a effective local interaction picture and have shown how this matches 
onto the usual low nuclear physics description of the process. 0] We demon- 
strate below how these two discussions can be merged — by determining the 
bound state energy via the presence of a pole in the scattering matrix — and 
thereby make contact with the traditional quantum mechanical discussion 
of hadronic atom energy shifts. |4j] In the next section we briefly review the 
usual quantum mechanics approach to the problem, while in section 3 we 
examine the same problem within the effective field theory procedure and 
demonstrate equivalence to the quantum mechanical results. In section 4 we 
discuss applications to pionium and to the ir~p atom. Finally, in section 5 
we present a short summary. 

2 Hadronic Atom Energy Shifts: Quantum 
Mechanical Approach 

The problem of calculating strong interaction energy shifts in hadronic atoms 
is an old and well-known one, and the methods by which to approach the 
subject are fairly standard. || In cases such as light pionic atoms or pionium 
the system is essentially nonrelativistic so that simple quantum mechanics 
can be employed. However, before considering such systems, we first examine 
the closely related problem of scattering of particles having the same charge, 
as considered by Kong and Ravndal. [0] 
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Consider then a system consisting of a pair of particles A, B both hav- 
ing charge +e and with reduced mass m r . Also, suppose that there exists 
no coupling to any other channel. First neglect the Coulombic interaction 
and consider only strong scattering. For simplicity, we represent the strong 
potential between the particles in terms of a simple square well of depth Vq 
and radius R — 



Considering, for simplicity, S-waves the wavefunction in the interior and ex- 
terior regions can be written as 



where jo, n-o are spherical harmonics and the interior, exterior wavenumbers 



are given by k = ^2m r E, K = ^2m r (E + V ) respectively. The connection 

between the two forms can be made by matching logarithmic derivatives, 
which yields the result 




(2) 




(3) 




k cot S = — — 1 + 



KRF(KR) 



1 



with 





x 



Making an effective range expansion 



k cot 5 — 

a 



+ . . . 



(5) 



we find an expression for the scattering length 



ao = R 1 — 



tan(K Q R) 

KnR 



where 





For later use we note that this can be written in the form 





The corresponding scattering amplitude is 



Now restore the Coulomb interaction in the exterior region. The analysis 
of the scattering proceeds as above but with the replacement of the exterior 
spherical Bessel functions by appropriate Coulomb wavefunctions Fq~, Gq 



jo(kr) 



n (kr) 



(9) 



whose explicit form can be found in reference || . For our purposes we require 
only the form of these functions in the limit kr « 1 — 



Kir) ^ «))(! + - + ...; 



Gt(r) 



kr«l 



+ 



C(r)+(k)) [kr 



a B 
1 



2 V+ (k) 



h{r 1+ {k)) + 2 1 



In 



2r 
a B 



+ 



Here 7 = 0.577215.. is the Euler constant, 

C\ x ) = 

exp(27ra) 



(10) 



(11) 



is the usual Coulombic enhancement factor, a B = l/m r a is the Bohr radius, 
fj+(k) = l/ka B , and 



h^ + (k)) = ReH(z7 1+ (k))=r 1 l(k)Y: 



\ n(n 2 + 7]\{k)) 



-\nr ]+ (k)- 1 (12) 



where H(x) is the analytic function 

H(x) = ijj(x) 



2x 



ln(x) 



(13) 



Equating interior and exterior logarithmic derivatives we find 
cos 5 Ff'(R) - sin 8 G£'(R) 



KF(KR) 



cos SqF^(R) - sin 5 G${R) 



kcot5 C\ V+ (k))^-^ 



k cot 6 C*(r l+ (k)) + i + ^ h(r, + (k)) - In ff + 2 7 - 1 



(14) 
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Since R « as Eq. O can be written in the form 



k cot 5 C 2 (r] + (k)) + — 
a B 



M^))-ln|| + 2 7 -l 



— (15) 
a 



The scattering length ac in the presence of the Coulomb interaction is con- 
ventionally defined as[[7| 

k cot 5 C 2 ( V+ (k)) + —h(rj + (k)) = -— + ... (16) 

ac 

so that we have the relation 

-- = (In -f + 1-27) (17) 

between the experimental scattering length — ac — and that which would exist 
in the absence of the Coulomb interaction — ao- 

As an aside we note that ao is not itself an observable since the Coulomb 
interaction cannot be turned off. However, in the case of the pp interaction 
isospin invariance requires a^ = a^ n so that one has the prediction 

-^=-^- aAWln SAb +1 - 27) < 18) 

While, strictly speaking, this is a model dependent result, Jackson and Blatt 
have shown by treating the interior Coulomb interaction perturbatively that 
a version of this result with 1 — > 0.824 is approximately valid for a wide range 
of strong interaction potentials || and the correction indicated in Eq. [18| is 
essential in restoring agreement between the widely discrepant — Oq™ = —18.8 
fm vs. a P c = —7.82 fm — values obtained experimentally. 

Returning to the problem at hand, the experimental scattering amplitude 
can then be written as 

f+(k) = e^>C\ V+ (k)) 



a iB h(i] + (k)) - ikC 2 (ri + (k)) 
e^C\ V+ (k)) 



where a = argr(l — irj + {k)) is the Coulomb phase. 

Analysis of the situation involving particles of opposite charge is similar, 
except that in this case in the absence of strong interaction effects there will 
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exist, of course, Coulomb bound states at momentum k n = in n = i/na B and 
energy E n = —K 2 l /2m r = —m r a 2 /2n 2 with n = 1, 2, 3, . . . In the presence 
of strong interactions between these particles, however, the energies will be 
shifted. One approach to the calculation of this shift is to examine the cor- 
responding scattering process A + B — > A + B. Then the existence of a 
bound state is indicated by the presence of a pole along the positive imag- 
inary axis — i.e. for n > under the analytic continuation k — > in. (If we 
(temporarily) neglect Coulomb effects and look only at strong scattering, this 
occurs when k = 1/ao in the case of Eq. ||.) 

However, in the case of oppositely charged particles A, B the analysis of 
the scattering amplitude must be in terms of appropriate Coulomb wavefunc- 
tions. If, as before, we include Coulomb effects only in the exterior region, 
then the appropriate forms of the wavefunctions for kr << 1 are given in 
ref.lEl as 



F "(r) k ^ C(r, + (k))(l 



+ 



a B 



kr«l 



C{r] + {k)) \ kr 



TV 2/ 1 

H(ir)+(k)) - i- coth(nr] + (k)) + 27 - 1 + lnz — 
2 ob 



+ . . . 

(20) 



Equating interior and exterior logarithmic deriviatives as before, we find 
cos SqFq'(R) — sin 5 Gq '(R) 



KF(KR) 



cos 5qF (R) — sin 5qG (R) 



-k C ot5 C\- V+ {k))^-^ 



k cot 6 C*(- V+ (k)) + i - %(h(ri + (k)) - ln|| + 2 7 - 1 



Thus we have 

k cot 5 C 2 (-r] + (k)) 



2 

a B 



%+(*))- la H + 27-1 



1 

a 



(21) 



(22) 



In this case the scattering length ac in the presence of the Coulomb interac- 
tion is defined as 



kcotSoC 2 ^-^)) h(rj+{k)) = 



a B 



(23) 
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so that we have the relation 



1 = -- + -(ln^ + l-2 7 ) (24) 



The corresponding scattering amplitude is then 

e- 2iao C 2 (-T] + (k)) 



fo(k) 



J- + JL h (rj + (k)) - ikC 2 (- V+ (k)) 



e- 2ia °C 2 (-ri + (k)) 



■± + ±[H(i7i+(k)) - Z7TCOth7T7 ?+ (fc))] 



(25) 



Under the continuation k — ► %k we have 

H{ir]+{k)) - in coth(nr] + (k)) -> /f (£) + vr cot vr£ (26) 
where £ = 1/kclb- The existence of a bound state is then signalled by 

+ — (#(0 + 7rcotO = (27) 

ac a# 

In the limit of no strong interaction — ac — > — we find £ n = l/K n dB = n 
and the usual Coulomb bound state energies 

* = -£ = -i£- " = 1 - 2 - 3 -- (28) 



while if ac* 7^ there exists a solution to Eq. |27 



Z = — *n + ^ (29) 

and a corresponding energy shift 

A£ n = -E n ^ + 0{^f (30) 
na% clb 

which is the conventional result. 0,[B| 

It is important to note here that Eq. ^ is written in a form that relates 
one experimental quantity — the energy shift AE n — to another — the scatter- 
ing length ac- Hence it is model-independent, even though, for clarity, we 
have employed a particular model in its derivation. This feature means that 
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it is an ideal case for an effective interaction approach, as will be shown in the 
next section. However, we first complete our quantum mechanical discussion. 

An alternative approach involves the use of bound state perturbation 
theory. |l| In this case the problem simplifies because of the feature that the 
range of the strong interaction — R — is much less than the Bohr radius — a#. 
Thus we may write 

AE n = <y\V\V>+Y / ^ + ••• 

~ |^(0)| 2 x lim < (f) f \V\^ +) > (31) 
Connection with the scattering length may be made via|| 

/W = -J<« (+) >^-« (32) 



so that for weak potentials we have 

m r f , q , r , s m r 4 
a ps 

in agreement with Eq. |^. Then since for Coulombic wavefunctions 



f d*rV(r) = -^7rR 3 V (33) 



l^(0)| 2 = ^r (34) 

we have 

AE n = #0 ) 2 + ... = -E n + ... 35 

m r na b 

as found via continuation of the scattering amplitude. (However, in this 
form it is not completely clear whether the relevant scattering length is the 
experimental quantity a>c or its Coulomb subtracted analog a .) 

In the real world, of course, this simple model is no longer valid, since re- 
alistic hadronic atoms, such as pionium or tt~p, are coupled to unbound 
systems — 7r°7r° or ir°n, — and must be treated as multi-channel prob- 
lems. Nevertheless the methods generalize straightforwardly from those given 
above. Specifically, in the absence of Coulomb interactions the scattering am- 
plitude for the 7T7T system can be given in the two-channel K-matrix formQ 

r^- A - 1 - 1 ( o 1 I ) < 36 > 

1 There exists also a coupling of tt + it~ via the anomaly to the ir°j channel, but this is 
p-wave and can be neglected for s-states such as considered here. 
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where ki,k 2 are the center of mass momenta in the 1 (neutral), 2 (charged) 
channels respectively and 

A =(T t) w 

\ a 21 a 22 ) 

is a real matrix given in terms of the coupled channel scattering lengths. The 
inverse of the matrix A is easily determined and the resultant form of the 
scattering amplitude is 

i 

/ 



1 + ik\a\\ + ik 2 a 22 — kik 2 detA 
( — an(l + ik 2 a 22 ) + ik 2 a 12 a 21 —a 12 
\ —a 2 i — a 22 (l + ikian) + ikiai 2 a 2 i 



(38) 



The existence of a bound state is, as before, indicated by the presence of a 
pole — i.e. by the condition 

1 + ikidu + ik 2 a 22 — kik 2 detA = (39) 

after approporiate analytic continuation. In the case of hadronic atoms, 
we must, of course, correctly include the Coulomb effects. From our single 
channel experience above, this is done via the prescription 

2 

ctij -> a£, -ik 2 -> — (H(£) + 7i cot vr£) (40) 
a B 



whereby the bound state condition — Eq. [39] — reads 

= 1 + -(H(0 + vr cot vrO (-ag + f^f °| ) (41) 

Thus we find 



AE n = -E n ^-(ag - ) + 0(^) 2 ( 12) 



as the coupled channel generalization of Eq. The real component of 

the energy shift is, of course, to lowest order identical to the single channel 
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result. What is new is that there has developed an imaginary component, 
corresponding to a decay width 



-2ImA£ n = -En—hlagl 2 + ... = -i^A;i|af 2 | 2 |*(0)| 2 + . . . (43) 
na B m Y> 



which is precisely what is expected from Fermi's golden rule 

r d^k-i k^ 



— 1\ < >j\ t | u-, -> | 2 

2n n^ml 1 ' 



(1) (2) 7T 

rrir mi- 



< 2 )— ^il*(0)r (44) 



or from second order perturbation theory. 



3 Hadronic Atom Energy Shifts: Effective 
Field Theory 

Identical results may be obtained from effective field theory and in many ways 
the derivation is clearer and more intuitive. ]10j First consider the situation 
that we have two particles A,B interacting only via a local strong interaction, 
so that the effective Lagrangian can be written as 

C = £ + - C tf A * A tf B * B + ... (45) 

The T-matrix is then given in terms of the multiple scattering series shown 
in Figure 1 

Tfi(k) = ~f{k) = Co + C 2 G (k) + C 3 Gl(k) + ...= i _ °° (46) 

where Go(k) is the amplitude for particles A, B to travel from zero separation 
to zero separation — i.e the propagator Dp(k; r* = 0, r = 0)— 

Go(fc) - ^oi (2vr)3 _*L _ + ie ~ J (2tt)3 P - s * + *e (4?) 

v ' 2m r 2m r 
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+ 




+ 



Figure 1: The multiple scattering series. 




+ 



Equivalently Tfi(k) satisfies a Lippman-Schwinger equation 

T fi {k) = C + C G (k)T fi (k). 



(48) 



whose solution is given in Eq. |46|. 

The function Go(k) is divergent and must be defined via some sort of 
regularization. There are a number of ways by which to do this. We shall 
herein use a cutoff regularization with k max = fi we have 



G (k) = -^(^+ik) 

Z7T 7T 



(49) 



Equivalently, one could subtract at an unphysical momentum point, as pro- 
posed by GegeliajlTl 



G (k) 



d 3 .< 



2m r 



+ 



2m r 



(2tt) 3 k 2 — s 2 + ie /i 2 + s 2 



-^(jM + ik) (50) 



which has been shown by Mehen and Stewart fT2|j to be equivalent to the 



PDS sceme of Kaplan, Savage and Wise.[T(J In any case, the would-be linear 
divergence is, of course, cancelled by introduction of a counterterm, which 
renormalizes Cq to Cq(/x). The scattering amplitude is then 



/(*) = -£ 



1 



G (k) 



2- 



C (m) " u ^/ m r C (M) 

Comparing with Eq. § we identify the scattering length as 

1 2tt 2/i 



^ -ik 



«o 



m r Co(fj,) 7r 



(51) 



(52) 



More interesting is the case where we restore the Coulomb interaction 
between the particles. The derivatives in Eq. ^5] then become covariant and 
the bubble sum is evaluated with static photon exchanges between each of 
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Figure 2: The Coulomb corrected bubble. 



the lines — each bubble is replaced by one involving a sum of zero, one, two, 
etc. Coulomb interactions, as shown in Figure 2. 

The net result in the case of same charge scattering is the replacement of 
the free propagator by its Coulomb analog 



G (k) - G+(*) = hm / 



d 3 s ipg (r')ipg* (f) 
(2vr)3-^--^+^ e 



d 3 s 2m r C 2 (r] + (s)) 
(2vr) 3 k 2 -s 2 + ie 



(53) 



where 

^+(f) = C{7 1+ {s))e lao e lS -\F 1 {-i7 1+ {s), l,isr - is- r) (54) 

is the outgoing Coulomb wavefunction for repulsive Coulomb scattering. ||13|| 
Also in the initial and final states the influence of static photon ex- 
changes must be included to all orders, which produces the factor 
C 2 {2nrj + {k)) exp(2io"o). Thus the repulsive Coulomb scattering amplitude 
becomes 

t+ru\- m r- CoC 2 (^+(fc))exp2iao 

tc {k) ~ ~ 2^ 1 - C Q G + C {k) (55) 

The momentum integration in Eq. can be performed as before using cutoff 
regularization, yielding 



Itx I 7r as 



H(iri + (k)) -hi^-( 



(56) 



where ( = ln27r — 7. (Equivalently, in the unphysical momentum subtraction 
scheme 

m 2 1 
G&k) = --l-(H(z V+ (k)) - H{ )) 

777 9 

~ --±<ji+ — [H(i V+ (k)) - In »a B -^(l)\) (57) 

Z7T 0>B 
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We have then 



C 2 (r] + (k))e 2tr7 ° 



2tt 



mrCo(jJ.) 



2/' 



2 



H(ir] + (k)) -ln^-C 



C 2 (v + (k))e 



2iao 



J_ 

ao 



2 



h( v+ (k)-\n^-(\ -tkC 2 ( V 4k)) 

Comparing with Eq. |19] we identify the Coulomb scattering length as 

1 _ 1 



(58) 



— ( ln ^T- + 



(59) 



which matches nicely with Eq. [T7| if a reasonable cutoff fi ~ ~ 1/i? is 
employed. The scattering amplitude then has the simple form 



C 2 ( V+ (k))e 



2iaQ 



(60) 



in agreement with Eq. |l9| . 

Now consider oppositely charged particles. In this case the analysis is 
parallel to that above, but there exist important new wrinkles in that the 
intermediate state sum in the Coulomb propagator must now include bound 

states 



Go(fc) - G~ c {k) 



lim 

r*,r— >0 



E 



fc 2 

nim 2m r 



i m r tt' ! 
~+~ 2n 2 



2m r 00 



E 



+ 



7raB ~^ n(n 2 + 77+ (A;)) V (27r) 3 /c 2 — s 2 + ie 



2m. 

d 3 s 2m r C 2 ( 



-»7(*)) 



where 



^j(f) = C(-r/ + (s))e- J<T0 e Jsr 1 F 1 (^ + (s), l,isr - is- r) 



(61) 



(62) 



is the outgoing Coulomb wavefunction for attractive Coulomb scattering and 



^ntmif) 



(2a 



B 



, (n-l-1)! 
2n[(n + £)\f_ 



e-^{2a B rYL 2 ^l l {2a B r)Y l 



(63) 



is the bound state wavefunction corresponding to quantum numbers nlm 
with Lj(x) being the associated Laguerre polynomial. Using the identity 



C 2 {-7] + {k)) = -C 2 (ri + (k)) + 2nr] + (k) cothnri+(k) 



(64) 
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we can write Eq. BT] as 



G~ c {k) 



-G%(h) + 2m r 



d 3 s 2nr] + (s) coth.7irj + (s) 



(27T) 



k 2 — s 2 + ie 



2 7T~bir> r ^ 

— E 



na B ~ n{n 2 + r]l(k)) 
m r (2/j, 2 
2% I 7T a# 



H(ii] + (k)) — in coth ni] + (k) — In 



2tt 



(65) 



where the integration is done via contour methods and the the contribution 
from the hyperbolic cotangent poles precisely cancels the bound state term. 
The resulting attractive Coulomb scattering amplitude is given by 



C 2 (- V+ (k))e 



-2icrQ 



2- 



2fi I 2 

m r Co(£t) tc ' a B 



H(irj+(k)) -incothnr] + (k) - In ft- - ( 



C 2 (-7] + (k))e- 2ia ° 



.i + JL 

a a s 



In ^-C] -i/cC 2 (— ?7 + (A;)) 
Identifying the Coulomb scattering length via 

1 1 2 aa B , 
= (In— — + C) 

<2<7 a a B ^ 



(66) 



this reduces to the simple form 



fc(k) 



C 2 (-7] + (k))e- 2iao 



i + ^[ H ^V + (k)) - iTTCOthTT^)] 



(67) 



(68) 



in agreement with Eq. In order to go to the bound state limit we can 
utilize the continuation 



H(ir] + (k)) — in cothnr] + (k) ^> H(£) + n cot 7r£ 
in which case we find the condition 





+ — (#(£)+7rcot7r£) 

a c a B 



(69) 



(70) 



in complete agreement with Eq. |27]. 
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In the real world — coupled channel — case the analysis is similar. We must 
generalize the effective Lagrangian to include off diagonal effects, so that Co 
becomes a matrix (Co)ij with i,j = 1,2 and the bubble sum must include 
two forms of intermediate states — charged and neutral. In the absence of 
Coulomb effects we have then a set of equations 

Tij = {Coh + J2( C ohG (k e )T ej (71) 

i 

whose solution is 

T= (l-C G (A;))- 1 Co (72) 

Explicitly, defining 
D = (1 - (C )iiGo(A;i))(l - (C ) 22 G (k 2 )) - {C^C^G^G^) (73) 
we find 

Tn = ((Co)n(l - (C ) 22 G (k 2 )) + (C Q ) 12 (C ) 21 G (k 2 ))/D 
T21 = (Cohi/D 
T 12 = (C ) 12 /D 

T 22 = ((Co) 2 2(l-(Co)iiGo(A;i)) + (Co)i2(Co)2iGo(A;i))/£> (74) 
Using the relation between the scattering amplitude and the T-matrix 



h = -M^T, (75, 

the connection between the parameters (C )jj and the corresponding scat- 
tering lengths dij is easily found via 

777(1) m (2) 



m (l) m (2) 
Ol2 = ^ —(C (ll)) 12 /J 



(1) (2) 

2 1 = — {C {n)) 2 i/J 

777(2) 777(1) 

022 = ^[(G)(/i))22-^detCb(/i)]/J (76) 
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with 



<rri(-0 777 777 (-0 777^) 

J=l- M^MCo(/i))n + ^(C (m))»] + /i 2 ^- ^- det C„(/0 (77) 

Z7T Z7T Z7T Z7T 

Then Eqs. [7B| are seen to be the same as the K-matrix forms Eg. |58} Inclusion 
of the Coulomb interactions is as before, involving the modifications 



T — > STS, where S 



1 

C(— rj + (k)) exp — ia 



2 

a ij ~^ a ?i an d — ik 2 — (H (£) + n cot ir£) , (78) 

a B 

and the resulting bound state condition is identical to Eq. [|l|. Thus the 
equivalence of the quantum mechanical and effective interaction methods is 
explicitly demonstrated. 

3.1 Effective Range Effects 

It is straightforward to go to higher order by inclusion of effective range 
effects. In the single channel quantum mechanical formulation this is accom- 
plished by the modification 

a c -> ac(l + ^a c r E k 2 + . . .) (79) 



Then Eq. [27] becomes 



a c 2£ 2 a| a B 

whose solution is 



' '' E + — (#(£) + Ti cot ttO ^ (80) 



AE n = —E n — — 



&B n \ a B a B J 



(81) 



so that effects are negligible. Similarly in the two channel case we alter Eq. 



via. 



Oij o»iU + ^ a ij( r E)ij(k 2 + fcj) + . . .) (82) 
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and the energy shift and width are modified appropriately 



ReA£ n = -E n ^- 
naB 



1 + a jn twin) ~-) + (^) 2 , 

Ob n \ as sa B 



T n = 2lmAE n = -E n ^k 1 \ag\ 2 (l + lag(r E ) 12 (k 2 1 -K 2 n )-(ag) 2 k 2 1 + ...) 

naB 2 

(83) 

Again for the real component the effective range effect is tiny and can gener- 
ally be neglected. However, in the case of the decay width there are two types 
of corrections which should be included — one due to the effective range cor- 
rection and a second due to multiple rescattering effects in the intermediate 
state. 

Using effective field theory the change is also directly obtained. In the 
single channel case, one begins by modifying the effective Lagrangian via 

£ -> C + ~C 2 (¥ A V* A ■ & B VV B + (Wjj^ • (W^b) (84) 

The multiple scattering series may summed exactly as before and the solution 
is found as(nj 

_ f C R (f,) + 2k 2 C 2 R (u))C 2 (-r 1+ (k))exp2 i a 
M J 1 1 - (C » + 2h*C?(ji))Ga(k) ^ 
where Cq (/z), C 2 (fi) are the renormalized quantities 

R _ C + LC\ amrfx „ R R _ C 2 - \ JC\ 

G " (1 _ C2 J)2 " ' °2 " {l _C 2j) 2 W>) 

and, in cutoff regularization, 

d 3 s „ 0/ , „ , /" <i 3 s 



7 = / #^ ( -^ +(s)) ' L = y (2^ s2c2( -^ +(s)) (87) 



Thus the scattering amplitude becomes 

C 2 {-r] + {k))e- 2ia ° 



fc(k) = ~T fi (k) 



2tt 2jx I 2 

"m r (C^(At)+2fc 2 C*| < (/i)) » + a B 



H(ir) + (k)) - in cothnr] + (k)) -ln^ff -£ 
C 2 (-77 + (A;)) e - 2 -° + ___ 
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Thus we identify 

l_ 2tt 2fi 2 87rC 2 /? (/i) 

and the bound state condition becomes Eq. ^(J, as expected. As pointed 
out in ref.||, the expression for effective range is not affected by Coulombic 
corrections and consequantly in the case of the NN interaction one expects 
equality for r^f and found experimentally. 

The coupled channel analysis is similar — the momentum dependent ef- 
fective coupling C 2 becomes a 2 x 2 matrix and the associated Lippman- 
Schwinger equation is algebraically somewhat more complex. However, the 
solution is straightforward and amounts in the end simply to a modification 
of the effective scattering length parameters defined in Eq. |76| by effec- 
tive range range effects as given in Eq. [8^, where (r^)^- are given in terms 
of the matrix elements C 2 ([i)ij. Thus the effective interaction bound state 
condition is identical to Eq. [18], as expected. 



4 Applications 

While the analysis given above is in some ways old, the applications are not 
and include two of the most interesting ongoing measurements in contempo- 
rary physics. One is the problem of pionium — 7r + 7r~ — whose existence has 



been claimed by a Russian group [|i~5f and which is now being sought in a 
major program at CERN. In this case then the first channel represents 7t 7t° 
while the second is 7r + 7r~. The reduced mass is 

m« = ^, m^ = ^f (90) 
and scattering lengths are given in terms of those with total isospin 0,2 via 

a n = ^(2ct2 + a ), a 12 = a 21 = ~^( a 2 - o ), a 22 = ^ (o 2 + 2a Q ) (91) 
where 

a ° = -*mf? a2= m^l (92) 

are the usual Weinberg values with F n ~ 92.4 MeV being the pion decay 
constant. [|16|] To lowest order we have a ground state energy shift 

Ait 

AisgVO = — (a 2 + 2a )|vl>(0)| 2 (93) 
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and a width 

rgVO = ^-J2m^m^a 2 - « ) 2 |^>(0)| 2 (94) 

which agree with the usual forms. M Corrections from higher order effects 
are found to be 



AE^J a B 



' !22 (2H(l)-3) 



r« i 

(2m n Am 7r )(-a 12 (r E ) 12 - a 2 u ) (95) 



9« 



(0) - " 7T,V 2 

Using the lowest order chiral symmetry predictions [| 



? uWu = 4fr ° ?1 = (96) 



we find a decay rate correction 



r(i) 
r (o) 

1 gs 



13% (97) 



in agreement with the result of Kong and Ravndal.|I|] Q In the case of the 
energy shift we find a negligible change 

AE^ 

£r = -8.2 x 1(T 4 (98) 

aeP 

(Strictly speaking, the charged channel phase shift here — a 22 — should be re- 
placed by its Coulomb corrected value 



---^ = - + -(m^f + l- 2 7 ) (09) 

a 22 a 22 a 22 a B *-H> 



However, this is only a small correction numerically.) The future detection 
of pionium should allow a relatively clean measurement of the irn scattering 
lengths. 

2 The dominant effect here is from the effective range, while the rescattering provides 
only a small correction. The form of our rescattering term differs from that of Kong and 
Ravndal, as we include only rescattering from physical — 7r°7r° — intermediate states while 
they include also that from (unphysical) charged states. 
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The other hadronic atom of current experimental interest is 7r~p wherein 
a PSI collaboration has already announced a measurement of the energy 
shift and decay rate and which plans to further improve these already pre- 
cise numbers. |L7| In this case we are dealing with a three channel system — 
n~p, Tr°n, 771. However, both the quantum mechanical calculation as well as 
the effective interaction analysis generalize and the bound state condition — 
Eq. |4l] — becomes (after some algebra) 

= 1 + —(#(£)+ 7T cot 7r£){-ag 
a B l 



ik\agag + ^3°32 a 23 ^1^3 a 23( a ll a 32 °31 a 12) + a 21 ( a 12 a 33 a 13 a 32 



+ 



C(„C„C „C„C\ 1 „C(„C„C „c„c 



1 + ifeian + iha-33 - hk^a^a^ - agag) 

(100) 



To lowest order then we have the results 

= -Ex— (hagag + hagag) (101) 

In the case of the energy shift, the result is as before, but in the case of the 
decay width we now have two contributions, one from the decay to the 7i°n 
channel and one from decay to 771. However, the tcN scattering lengths can 
then be obtained once the radiative channel is subtracted off via 

r.o n = (102) 



where P = T-n/T^On = 1.546 ± 0.009 is the Panofsky ratio. ||T8 



On the theoretical side the 7tN scattering lengths can be written in terms 
of those with total isospin |, | via 

1 v^2 1 

a n = + 2a 3 ), aw = a 21 = — (a 3 - at), a 22 = -(2% + a 3 ) (103) 

and the lowest order chiral Lagrangian yields values for these quantities 

777™ 1 

«■ = = -S%TT^ (104) 
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and higher order chiral corrections have also been calculated. [HJ In the case 
of the radiative channel, the scattering length is also predicted by chiral 
symmetry in terms of the Kroll- Ruder man term[f20f 



where g& ~ 1.25 is the axial coupling in neutron beta decay. 

There exist higher order corrections to these lowest order predictions from 
both rescattering effects — cf. Eq. |100| — and from inclusion of effective range 
corrections. However, unlike pionium the latter are not given by chiral sym- 
metry and hence are model dependent. Nevertheless when known corrections 
are included present experimental results extracted from n~p experiments 

6i = Ua 3 - ai ) ~ 0.096 ± 0.007m; 1 , b ~ h ai + 2a 3 ) = 0.0105 ± 0.007m; 1 

(106) 

are in satisfactory agreement with the chiral symmetry predict ionQ 

0.096m; 1 < 6i < 0.088m; 1 (107) 

although there is at present a small discrepancy with a parallel measurement 
of the shift in deuterium. [|H]] However, our purpose here is not to discuss 



experimental interpretations, rather only to point out the connection between 
conventional and effective interaction methods. 



5 Conclusions 

We have above analyzed the problem of hadronic atom energy shifts due to 
strong interaction effects both via a traditional quantum mechanical discus- 
sion and via a calculation in the effective interaction picture wherein the low 
energy hadronic interactions are written in terms of an effective local poten- 
tial and have demonstrated the complete equivalence between the two pro- 
cedures. Applications have been made to the systems pionium — 7r + 7r~ — and 
tt~P currently being studied experimentally at CERN and PSI respectively. 
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